Sub-micrometer epitaxial Josephson junctions for 
quantum circuits 

Jeffrey S. Kline 1 , Michael R. Vissers 1 , Fabio C. S. da Silva 1 , 
David S. Wisbeym Martin Weides 1 , Terence J. Weir 2 , 
Benjamin Tureki§|, Danielle A. Braje 2 , William D. Oliver 2 , 
Yoni Shalibo 3 , Nadav Katz 3 , Blake R. Johnson 4 , Thomas A. 
Ohki 4 , and David P. Pappas 1 

1 National Institute of Standards and Technology, Boulder, CO 80305, USA 
2 MIT Lincoln Laboratory, 244 Wood Street, Lexington, MA 02420, USA 
3 Racali Institute of Physics, The Hebrew University of Jerusalem, Jerusalem 91904, 
Israel 

"Raytheon BBN Technologies, Cambridge, MA 02138, USA 
E-mail: j ef f rey . kline@NIST . gov, david . pappasQNIST . gov 

Abstract. We present a fabrication scheme and testing results for epitaxial sub- 
micromcter Josephson junctions. The junctions are made using a high-tempcrature 
(1170 K) "via process" yielding junctions as small as 0.8 ^m in diameter by use 
of optical lithography. Sapphire (AI2O3) tunnel-barriers arc grown on an epitaxial 
Rc/Ti multilayer base-electrode. We have fabricated devices with both Re and Al top 
electrodes. While room-temperature (295 K) resistance versus area data are favorable 
for both types of top electrodes, the low-temperature (50 mK) data show that junctions 
with the Al top electrode have a much higher subgap resistance. The microwave loss 
properties of the junctions have been measured by use of superconducting Josephson 
junction qubits. The results show that high subgap resistance correlates to improved 
qubit performance. 
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1. Introduction 

Josephson junction superconducting devices are promising candidates for qubits in 
quantum information circuits [I]. The tunnel-barriers in these junctions are typically 
amorphous A10 x made by room-temperature (T = 295 K) oxidation of thin films 
of aluminum. Qubit spectroscopy reveals coupling to stochastically distributed two- 
level systems (TLS) in the tunnel-barrier [2]. These TLS are observed as avoided level 
crossings (i.e., splittings) in the qubit spectroscopy. For the amorphous A10 x tunnel- 
barrier, the density of TLS splittings is measured to be 0.5 (/im 2 GHz) -1 [31 H]. While 
the physical origin of TLS is still under debate, it is clear that their interaction with the 
qubit is detrimental because they can absorb energy and decohere the qubit state. These 
TLS have a random distribution in frequency space and coupling strength. Unless some 
strategy for reducing the number of TLS is used, it is highly likely that TLS splittings 
will appear close to the desired operation frequency when circuits with multiple qubits 
are constructed. 

To date, there have been two strategies to reduce the number of TLS in qubit 
junctions. The first is to reduce the junction area as much as possible. Sub- 
micrometer Josephson junctions made by use of electron-beam lithography and Al 
shadow-evaporation [5] are highly successful in charge qubits, transmons, flux qubits, 
and low- impedance flux qubits [6j [7J El E] • However, the absence of metal cross-overs in 
electron beam-defined circuits limits the available circuit designs (e.g., no gradiometric 
flux coils). Step-edge junctions fabricated by use of optical lithography for phase qubits 
can be made as small as 1 fim 2 [101 [TTj . While cross-overs are part of step-edge 
technology, multiple qubit circuits will still suffer from the residual stochastic TLS 
splitting distribution due to the high density of TLS splittings. 

The second strategy is to reduce TLS density by use of epitaxial materials. 
While this typically involves high-temperature processing [12], it has yielded improved 
performance. Oh et al. observed an ~80% reduction in the density of TLS in a large- 
area (70 /iin 2 ) phase qubit with a crystalline AI2O3 tunnel-barrier when compared to 
amorphous A10 x [13]. This technology uses optical lithography, and cross-overs are 
easily made. 

Our goal in this work is to combine these two strategies by developing a process 
to reduce the size of epitaxial junctions for high- coherence qubits. In addition, we 
evaluate the efficacy of replacing the Al top electrode with Re. This is motivated by the 
discussion in [13], where it was hypothesized that the residual TLS may originate at the 
A1 2 03-A1 interface. To test this hypothesis, we studied Re top-electrode junctions and 
qubits, and compared them to devices with Al top-electrodes. 

In order to reduce the junction size, we first tried a standard trilayer process [T2] 
but the photoresist pillar washed away in the developer rinse for junction sizes smaller 
than ~2 /jm. While a trilayer process for sub-micrometer junctions does exist [14], 
it requires chemical-mechanical planarization and this is not available in our facility. 
Instead, we developed a high-temperature "via process" similar to a scheme used for 
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masking GaN nanowire growth [TB]. As discussed below, the epitaxial base-electrode 
is grown on a high-quality substrate, an insulator with a via to the base-electrode is 
defined, and the epitaxial tunnel-barrier is then grown in the via after heating and 
recrystallizing the surface. We have measured both the room-temperature (T = 295 K) 
and low-temperature (T < 100 mK) properties of these single-junction devices, and we 
have also fabricated qubit devices and measured their performance. 

2. Substrate preparation and base-electrode growth 

All of our devices were fabricated on single-crystal Al 2 03(0001) sapphire wafers. The 
wafers are 76.2 mm in diameter and 0.43 mm thick. As received from the manufacturer, 
the surface of the substrate exhibits no lateral crystalline order when imaged by atomic 
force microscopy (AFM), i.e., it is amorphous. To improve the crystalline order, we 
heat the substrate in a tube furnace to 1370 K for 20 hours in a 14:1 nitrogen-to-oxygen 
gas mixture at atmospheric pressure [161 02] • After the furnace treatment, we observe 
atomic step terraces and lateral order. In addition, we find a correlation between the 
miscut angle (angle between dicing saw cut and the (0001) crystal plane) and terrace 
size measured by AFM: even a small, 0.3° miscut limits the terraces to ~100 nm wide, 
while a nominal 0.0° miscut yields ~390 nm wide terraces. However, we found that the 
surface morphology of the Re base-electrode is independent of furnace treatment and 
miscut angle. 

At first we used a 165 nm thick rhenium film for the base-electrode, deposited 
by use of ultra-high vacuum (UHV) DC sputtering [18]. The substrate temperature is 
held at 1170 K, the deposition rate is 3 nm/min, and the argon sputter gas pressure is 
0.7 Pa. For this and all subsequent layers, we rotate the substrate during deposition 
to improve film thickness uniformity. The magnetron sputter guns are mounted in a 
sputter-up configuration at 35° off normal and 15 cm from the substrate. Using this 
technique, we obtain crystalline rhenium films. We find that these film are characterized 
by ~100 nm diameter hexagonal islands with ~15 nm height, as shown in figure [T](a). 
The root mean square (rms) roughness of the Re films is 3.2 nm, and is indicative 
of step-bunching and limited mobility of the Re during deposition. For comparison, a 
polycrystalline or crystalline Nb base-electrode suitable for high-quality Nb-Al/AlO^-Nb 
Josephson junctions has roughness < 0.5 nm [Tj)J 120] . 

In order to obtain a smoother Re surface for subsequent growth of the barrier, we 
found that it is possible to reduce the rms roughness of the base-electrode film while 
maintaining crystallinity by using a Re/Ti multilayer. In this process, we deposit a 
10 nm Re layer and then cap it with 1.5 nm of Ti. Both films are UHV sputtered at 
1170 K. Titanium has a lower surface free energy (1.9 J/m 2 ) than rhenium (2.2 J/m 2 ) 
and acts as a wetting layer, resulting in a significantly smoother surface, as shown in 
figures [Q (a) and (c). By repeating the Re/Ti unit cell structure twelve times and then 
capping with a 10 nm Re top layer (i.e., (Re/Ti) 12 Re), we obtain base-electrode films 
150 nm thick with an rms roughness of only 0.6 nm. Line scans from the respective 
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AFM images are shown in figures [T] (b) and (d), illustrating that the multilayer film is 
much smoother, with fewer vertical edges, than the pure Re film. Reflection high energy 
electron diffraction (RHEED) patterns (not shown) from pure Re and (Re/Ti)i2Re films 
are indistiguishable, indicating that base-electrode crystallinity is not degraded by using 
the Re/Ti multilayer. This base-electrode surface is much more favorable for tunnel- 
barrier growth with a sharp metal-insulator interface [191 EU 122] ■ 

3. Tunnel-barrier and top-electrode growth 

Once the base-electrode is grown, we proceed to define the tunnel junction and top- 
electrode by use of the via process. This process is illustrated in figure El The tunnel- 
barrier and top-electrode are deposited after the insulator and via are defined. To 
accomplish this, we remove the wafer (with the epitaxial base-electrode already grown) 
from the UHV sputter tool. The first step of the process is shown in figure E^a), where 
the (Re/Ti) 12 Re base-electrode is patterned by use of a 500 V SF 6 RIE at 2 Pa, etching 
all the way down to the substrate. The base-electrode is then covered with a 220 nm 
thick cross-over insulator, either SiO x or SiN^, by use of plasma enhanced chemical 
vapor deposition at T = 295 K. 

Vias are then etched in the insulator by use of a 240 V RIE, shown in figure |2^b). 
We use CHF 3 + 2 at 13 Pa for the SiO x insulator and CF 4 + 2 at 2 Pa for SiN^. 
This step defines the tunnel junction size and shape. This RIE has a 3:1 (5:1) selectivity 
in etch rate between SiO x (SiNj.) and Re, allowing us to stop fairly effectively when we 
reach the top of the crystalline Re base-electrode. Minimizing the over-etch into the 
base-electrode is critical, as it could create vertical walls around the edge of the via and 
it also could etch down to the Ti wetting layer. Tunnel-barrier coverage on these vertical 
sidewalls would be problematic and prone to pinholes and uneven coverage. We use a 
laser interferometer endpoint-detection scheme to minimize the over-etch (typically < 10 
% of the total etch time). We estimate the maximum amount of Re removed by the 
via over-etch and RF-clean (see below) is ~5 nm. This leaves ~5 nm of Re covering 
the nearest Ti wetting layer, so that our tunnel barriers are grown on Re, not Ti. This 
assertion is also confirmed by RHEED patterns: Re and Ti have distinct patterns. 

After the via etch, the base-electrode in the bottom of the via has been amorphized 
due to the over-etch portion of the via etch. This surface needs to be cleaned and 
recrystallized before a tunnel junction can be grown on it. We do this by loading 
the wafer back into the UHV sputter tool, performing an argon RF-clean to remove 
~2 nm of material, and then heating the wafer to 1170 K for 1 hour as shown in 
figure El^c). RHEED image shown in the inset of[2](c) and ex-situ AFM images (not 
shown) indicate that the Re surface is clean and re-crystallized. The epitaxial AI2O3 
tunnel-barrier is then grown in situ by use of UHV RF magnetron-sputtering from a 
sintered AI2O3 sputter target [231 El]- For this deposition, the substrate temperature 
is held at 1170 K, the deposition rate is 0.9 nm/min, and the sputter gas pressure 
is 0.7 Pa argon with 5 mPa oxygen. The oxygen gas is necessary to prevent oxygen 
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loss from the aluminum oxide at high-temperature and to obtain fully stoichiometric 
AI2O3. The thickness of the tunnel-barrier is monitored in situ by use of spectroscopic 
ellipsometry. We grow crystalline aluminum oxide films 1.8 ± 0.2 nm thick as tunnel- 
barriers. We find that they are conformal to the Re base-electrode, as evaluated by 
comparing AFM images and finding them to be indistinguishable from those of the 
base-electrode. In addition, we note here that other barrier growth conditions were 
explored, for example, growing at 1170 K without oxygen and also growth at T = 295 
K followed by an 1170 K anneal in oxygen. These resulted in a low resistance x area 
(RA)-product < 400 f2/im 2 for tunnel-barriers up to 9 nm in thickness, and the RA- 
product was independent of barrier thickness. From this, we conclude that tunnelling 
was not the dominant transport mechanism for barriers grown without oxygen or at 
T = 295 K with an anneal. According to AFM and electrical-isolation measurements of 
metallic cross-overs, the SiO x . and SiN x cross-over insulators are stable (i.e., no flowing 
of insulator material) and isolate well even after the 1170 K processing. 

The top-electrode (either Re or Al) is then deposited in situ by use of UHV DC 
magnetron-sputtering after the wafer is cooled to room-temperature (T = 295 K) in 
a 5 mPa oxygen background. The Al is deposited at a rate of 3 nm/min, and the 
argon sputter gas pressure is 0.7 Pa. For the Re top-electrode, we use xenon sputter 
gas instead of argon to avoid the creation of energetic neutral sputter-gas atoms, which 
act as an unintentional mill of the tunnel-barrier during the first few atomic layers of 
top-electrode deposition. This is because energetic neutrals are created when there is 
a large mismatch in atomic mass between the sputter gas and the target material [25J. 
The use of xenon instead of argon for Re sputtering reduces the fractional energy of 
neutrals from 0.42 to 0.03. According to RHEED and AFM (not shown), the Al and Re 
top-electrodes exhibit moderately textured in-plane crystalline order, but small ~ 30 
nm grain-size due to the low 295 K deposition temperature. In the final step, figure 
[2(d), the Re top-electrode is patterned by use of a 500 V SFg RIE at 2 Pa. If the 
top-electrode is Al, we use a 200 V argon ion mill at 0.4 Pa (oriented 20 degrees from 
substrate normal with sample rotation). 

4. Electrical Characterization 

We measure the room-temperature (T = 295 K) resistance of octagonal test junctions 
ranging in designed minimal diameter d from 0.5 /im to 15 /im (area: 0.2 fim 2 to 186 
/mi 2 ). This provides three important pieces of information, including the process bias 
do. First, for medium and large junction sizes (d » d ), the RA-product should be 
flat when plotted versus designed area if there are no spurious transport channels at 
the perimeter of the junction. Second, by plotting RA-product versus electrical area 
0.827(d — d ) 2 and adjusting d so that we obtain a flat RA-product curve for small 
junction sizes (d ~ do), we extract do- This gives us information concerning how the 
actual size of the junction differs from the designed size. Third, if the superconducting 
gaps of the top and base-electrodes are known, the critical current density for the 
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junctions in the superconducting regime can be calculated [2E]. This gives us feedback 
to adjust the RA-product by changing the tunnel-barrier deposition time for subsequent 
wafers. 

Figure [3] shows a plot of RA-product versus electrical area for junctions with Re 
and Al top-electrodes. For both types, the curve is flat for medium and large junctions, 
so we expect no significant perimeter transport. Both types of junctions have a process 
bias of -0.3 /im, meaning that the junctions are 0.3 /im larger in diameter than designed. 
This agrees well with the SEM image in figure 12(d), where a junction that was designed 
as 0.5 /im was measured to be 0.8 /im (area = 0.5 /im 2 ). In order to account for the 
observed spread in RA-product, we designed qubit circuits with various sized junctions, 
as described in [4]. Based on measurements at T = 295 K, both types of top-electrodes 
appear favorable for use as Josephson junctions. 

Low-temperature (T ~ 50 mK) measurements were then conducted for both the 
Re and Al top-electrode devices in an adiabatic demagnetization refrigerator using a 
commercial data-acquisition card and preamplifier. Figure H] shows IV curves for two 
junctions of nominally equal area (~60 /im 2 ) and RA product (~2000 Qfim 2 ). While the 
normal-state resistances, i.e., the inverse slope of the curves above the superconducting 
gaps, are nearly the same, a dramatic difference in the subgap structure is observed. 
For the Re top-electrode junction, we see low subgap resistance R sg = 226 Q, indicating 
some transport mechanism other than pure tunnelling. The subgap resistance is only five 
times higher than the normal-state resistance. The Al top-electrode junction shows a 
sharp corner, high subgap resistance and a re-trapping current that is limited by system 
noise, indicative of a high-quality junction [21]. We measured tunnel junctions ranging 
in size from 0.5 — 186 /im 2 from five wafers with Re top-electrodes and six wafers with 
Al top-electrodes: all measurements exhibit the same qualitative behavior where the Re 
top-electrode junctions have low subgap resistance and the Al top-electrode junctions 
have high subgap resistance. We conclude that junctions made using Re top-electrodes 
have inherently poor subgap properties. 

We measured the superconducting critical temperatures of the electrodes: 1.1 K 
(Al), 2.5 K (Re) and 2.4 K ((Re/Ti) 12 Re multilayer), corresponding to superconducting 
gaps A of 0.17 meV (Al), 0.38 meV (Re) and 0.36 meV ((Re/Ti)i 2 Re), using BCS theory 
[28]. The measured values Ai + A 2 of 0.75 meV for (Re/Ti)i 2 Re-Al 2 3 -Re and 0.45 
meV for (Re/Ti) 12 Re-Al 2 03-Al from the IV curves in figure H] are in good agreement 
with theory. 

We also measured superconducting qubits made using the via process with both Re 
and Al top-electrodes. We first describe a flux-biased phase qubit with Re top-electrode. 
The circuit design is similar to [4] with qubit state measurement performed using a DC 
SQUID. For a device with a 4 /im 2 qubit junction with capacitance ~ 200 fF, critical 
current = 2 //A, 700 fF shunt (Re/Ti)i 2 Re interdigitated capacitor, loop inductance 
L = 720 pH, and 1 fF SiO x . cross-over insulator, we measured an energy relaxation time 
Ti = 15 ns, as shown in figure [5(a). We hypothesize that T\ is limited by the relatively 
low subgap resistance of the qubit junction with Re top-electrode; the classical RC 
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decay time for the qubit is r = CR sg ~ 2 ns, where C = 900 fF is the total qubit 
capacitance. We measured two phase qubits from two wafers and both yielded similar 
results. We were unable to detect TLS splittings in the spectroscopy data due to the 
broad linewidth caused by the short relaxation time of these qubits. 

A qubit fabricated with an Al top-electrode showed a much longer T\ time of 500 
ns, as shown in figure [3£b). These data were taken from a transmission-line shunted 
plasma oscillation qubit (transmon) with dispersive qubit state readout using a half- 
wave resonator [29|. The total qubit capacitance is given by two 1 fim 2 junctions with 
~ 100 fF capacitance (critical current = 0.1 /xA), 60 fF shunt (Re/Ti) 12 Re interdigitated 
capacitor, and 1 fF cross-over SiN x insulator. The half-wave resonator frequency is 8.3 
GHz and the T\ measurement was performed at the 7.3 GHz flux "sweet spot". We 
measured two transmon qubits from one wafer and both yielded similar results. Based 
on the qubit-resonator coupling of g/lix = 85 MHz, qub it-resonator detuning A/2n = 1 
GHz, and resonator photon loss rate k/2tt = 0.8 MHz: the Purcell effect limit on T\ is 
(A/g 2 )/hi = 27 /is, so our devices are not limited by the Purcell effect. We observed 
three TLS splittings in the spectroscopy measurement over a 0.5 GHz range (not shown), 
with maximum splitting size = 7 MHz. 

Table [1] shows an analysis of the loss in each element of the transmon circuit: 
Josephson junction, interdigitated capacitor, and SiN x insulator. The participation ratio 
of each element is given by pi = Ci/C to t, where Cj is the capacitance of element i and C to t 
is the total qubit capacitance. The contributed loss is given by the microwave dielectric 
loss tangent (tan 5) times P«. Here we consider only the low-power loss tangent, i.e., 
the loss tangent measured when the TLS are unsaturated by the applied electric field 
[3] . We use independently measured values of tan S for the interdigitated capacitor and 
the SiNa, insulator. We use the measured T\ = 500 ns to calculate the total loss tangent 
of the transmon as 4.3 x 10~ 5 through 7\ = l/(27r/ r tan S), where f r is the 7.3 GHz 
resonance frequency. We find that the performance of the qubit is limited primarily by 
loss in the Josephson junction and the interdigitated capacitor. Other loss mechanisms, 
such as non-equilibrium quasiparticles, are not considered in this analysis. 

Table 1. Transmon loss analysis. The capacitance of element i is Ci, the participation 
ratio of element i is Pi, the loss tangent is tan 8, and the contributed loss is given by 
tan S x Pi. 



Element 


a (fF) 


P* (%) 


tan (5 


Contributed loss 


Junction 


100 


62.1 


3.5 x 10~ 5 


2.2 x 10~ 5 


IDC 


60 


37.3 


4.0 x 10~ 5 


1.5 x 10~ 5 


SiNz 


1 


0.6 


1.0 x 10" 3 


6.2 x 10~ 6 








Total loss 


4.3 x 10~ 5 
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5. Conclusions 

We have presented a recipe for the fabrication of sub-micrometer epitaxial Josephson 
junctions with Al 2 03 tunnel-barriers. The substrate crystallinity has been improved 
by a furnace anneal, and the base-electrode has been smoothed through the use of a 
(Re/Ti)i 2 Re multilayer base-electrode. The epitaxial Al 2 03 tunnel-barrier is deposited 
at the bottom of a via in either SiO x or SiN x . The top-electrodes are made from either 
Re or Al. 

We find that Josephson junctions fabricated using the via process with Re top- 
electrodes have low subgap resistance and phase qubit energy relaxation time T\ — 15 
ns. This energy relaxation time is much smaller than the T\ = 500 ns measured on a 
large area (49 /mi 2 ) epitaxial Re-Al 2 03-Al phase qubit [I] fabricated using a trilayer 
process in the same laboratory as the devices studied in this work and also the best 
amorphous-barrier phase qubit with 7\ = 600 ns [30]. We find that our Al top-electrode 
devices have a high junction subgap resistance and transmon qubit energy relaxation 
time T\ = 500 ns. We note that the best amorphous-barrier transmon, with junction 
area ~ 0.1 /im 2 , has T\ = 2000 ns for operation at / ~ 6 GHz [31] . 
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Figure 1. Comparison of base-electrode film morphology. Rhenium base-electrode 1 65 
nm thick, sputter-deposited at f f70 K: (a) AFM image (lxf fim 2 ) with rms roughness 
= 3.2 nm and (b) AFM linescan. (Re/Ti)i2Re multilayer base-electrode film 150 nm 
thick, sputter-deposited at 1170 K: (c) AFM image (lxl /im 2 ) with rms roughness = 
0.6 nm and (d) AFM line scan. Z-scale is 25 nm for both images. 
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Figure 2. Fabrication schematic of epitaxialJoscphson junctions, (a) the (Rc/Ti)i2Rc 
base electrode is patterned and encapsulated with insulator, (b) The via is etched into 
the insulator, (c) RF-clcan, 1170 K outgas, tunnel-barrier growth, and top-electrode 
deposition are all performed in UHV without breaking vacuum. The inset shows a 
crystalline RHEED pattern of a Re base-electrode after 1170 K anneal, (d) The top- 
electrode is patterned. Inset shows an SEM image of a junction with 0.5 ^tm designed 
diameter (measured diameter = 0.8 /im). Additional wiring layers are needed for 
gradiometric devices. 
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Figure 3. Room-temperature (T = 295 K) resistance measurements of tunnel 
junctions with Re (black ■ curve) and Al (red • curve) top-electrodes. Resistance 
measured at 100 nA bias current. 
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Figure 4. Low-temperature (T ~ 50 mK) I-V curves for the Josephson junctions 
with Re (black ■ curve) and Al (red • curve) top-electrodes. The superconducting 
branch is intentionally suppressed by a magnetic field oriented in the plane of the 
junction. We measure Ai + A2 = 0.75 meV for the junction with Re top-electrode, 
and Ai + A2 = 0.45 meV for the junction with Al top-electrode. 
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